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Abstract 

The problem of quantizing the radiation field inside a nonlinear dielectric is 
studied. Based on the quantization of radiation in a linear dielectric which 
includes absorption and dispersion, we extend the theory in order to treat also 
nonlinear optical processes. We derive propagation equations in space and 
time for the quantized radiation field including the effects of linear absorption 
and dispersion as well as nonlinear optical effects. As a special case we derive 
the propagation equation of a narrow-frequency band light pulse in a Kerr 
medium. 
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I. INTRODUCTION 



Non-classical properties of quantum light pulses propagating through nonlinear optical 
media are of increasing interest. E.g. quantum solitons in fibers may in future play an 
important role in communication technology ffl. But the proper consideration of quantum 
effects as well as changes of quantum properties of light pulses propagating through nonlinear 
optical media requires an adequate description of the quantized radiation field in such media. 
In addition to the nonlinear optical properties of the dielectric under consideration, also 
effects of linear absorption and dispersion must be included. As it is well known the existence 
of solitons in fibers with a Kerr-like nonlinearity requires also the dispersion of the group 
velocity of the wave [@,[3| . On the other hand non-classical phenomena such as squeezing are 
very sensitive to degradation effects as for instance absorption of radiation. Therefore it is 
necessary, to develop a theory of the quantized radiation field in dielectrics which considers 
both nonlinear optical effects and dispersion as well as absorption effects. The quantization 
of the radiation field in dielectrics has been performed in the past in different ways. Some 
of this treatments are presented in p]-2"T||. The quantization of the electromagnetic field in 
linear, lossless, and dispersionless dielectric media has been developed in |5|-[10|. The case 
of lossless and dispersionless nonlinear dielectrics has been treated by Hillery and Mlodinov 
TT] . Squeezing of solitons in dispersive nonlinear dielectrics was successfully explained in 



in 



|12|1 . A canonical quantization schema for a general dispersive nonlinear dielectric has been 



presented by Drummond |L5| . The methods of the phase-space formulation of the equations 
for quantum optical pulse propagation in nonlinear and dispersive single-mode fibers are 
discussed in . A very interesting approach to the quantization of the radiation field in a 



dielectric composed of two- level atoms has been recently proposed in pi|. 

Huttner and Barnett [16|] presented a canonical quantization scheme for the radiation 
field in a linear dielectric, employing the Hopfield model of the dielectric [|J. It is worth not- 
ing that in this approach both the dispersion and the absorption by the medium are taken 
into account in a quantum-mechanically consistent way. A phenomenological generalization 
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of this theory was given in fl8-20|,23|, and the theory has been applied to the propagation 



of quantum pulses in . In this paper we develop a theory of the quantized radiation field 
in dielectrics which includes both nonlinear optical effects and dispersion as well as absorp- 
tion effects. We employ the quantization scheme developed by Huttner and Barnett ||16| . 
This quantization schema is based on a microscopic model of Hopfield type ||, where the 
dielectric medium is represented by a polarization field X (x). We extend this quantization 
scheme by introducing nonlinear polarization field terms into the Hamiltonian. In reality 
the material system is not a linear one, and these additional terms describe the intrinsic 
nonlinear properties of the dielectric which can give rise to nonlinear optical effects. In 
Sec.[H] we review the main results of the quantization scheme for a linear dielectric described 



in [|16|,£(J . In Sec.|HI| the theory is extended to nonlinear dielectrics, and the quantization 
schema in this case is described, including linear dispersion and absorption. A propagation 
equation for a narrow-frequency band quantum pulse is derived in Sec.[TV| and the special 
case of a Kerr-like nonlinearity is considered in more detail. A summary is given in Sec.|V|. 



II. LINEAR MEDIUM 

In this section we review some results of the quantization of the electromagnetic field in 
a dispersive and absorptive linear dielectric. We employ the quantization scheme developed 
by Huttner and Barnett [TJ| and extended in ||20|| . This quantization schema is based on 
a microscopic model of Hopfield type [|J where the dielectric medium is represented by a 
polarization field X (x) , which is coupled to a reservoir composed of a continuum of harmonic 
oscillators Y (x, u) to allow for absorption (0 < uo < oo, the continuum of bath modes at 
a fixed space point x is labeled by the parameter a;). For convenience we consider only a 
one-dimensional model of the dielectric, the propagation direction of the radiation field is 
the x-direction and the vector potential A (x) is linear polarized in the z-direction. The 
electromagnetic field is coupled to the polarization field in the minimal coupling version, 
and the resulting Hamiltonian can be diagonalized introducing bosonic operators / (x, uj) 
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and P (x, uj), which describe polariton-like excitations inside the dielectric Hl6,20|. After 
diagonalization the Hamiltonian of the transverse electromagnetic field interacting with the 
linear dielectric can be written as 



H = hH L = h / dxH L (x 



(1) 



where 



H L (x) = / duo uj P (x, to) f (x, to) . 



(2) 

The creation and destruction operators p (x, uj) and / (x, uj) satisfy the familiar bosonic 
equal time commutation relations 



/ (x, uj) , p (x', uj') =5{x-x') 8{u- uj') 



(3) 



f{x,u)J{x',uj') 



0. 



(4) 



The operators of the electromagnetic field, e.g. the vector potential A (x) and the electric 
field strength E (x) as well as the polarization field operator X (x) and the bath operators 
Y (x, u) may be written in terms of the bosonic operators p (x, uj) and / (x, u) as follows 



A (x) = (x) + (x) = ( duo \ A (x, to) + i f (x, v) 



(5) 



E (x) = E (+) (x) + E ( - ] (x) = j duo[E (x, uj) + & (x, uj) 





(6) 



where 



X (x) = X {+) (x) + (x) = j duo [X (x, uj) + P (x, uj) 





A (x, uj) = / dx' Ga (x, x', uj) f (x', uj) 



(7) 



(8) 



E (x, uj) = iuoA (x, uj) 



(9) 



—X (x, oj) = (e (oj) - 1) E (x, oj) - 2iacJe~{u) f (x, oj) . (10) 
Co 



The Green function Ga (x, x', oj) is given by 

Ga (x, x', oj) = —ia, 



gi M ik(u)\x-x'\ nn 

V e (w) ' l j 



where £ (a;) = e r (oj) + ie\ (oj) denotes the complex permittivity of the dielectric. The complex 
wave number k (u) is related to the complex permittivity e (u) and the complex index of 
refraction n (uS) by 



k (u) = — J e (oj) = — n (oj) = — (n r (oj) + n\ (oj)) = k T + ik{. (12) 
c v c c 

A normalization constant a has been introduced for convenience 



tt = V4^4' (13) 

where eo is vacuum permittivity constant, and A is the normalization area perpendicular to 
the x direction. The coefficient p denotes the coupling constant between the polarization field 
X (x) and the electromagnetic field. It is worth noting that the polarization field operator 
X (x, oj) may be expressed as a sum of the electric field operator E (x, oj) multiplied by 
(e (oj) — 1) and the bosonic operator / (x, oj) multiplied by the square root of the imaginary 
part of the complex permittivity, which describes absorption in the dielectric. Explicit 



expression for the complex permittivity e (oj) can be found in [JT6], but these are not of 
relevance for our task, and any phenomenologically introduced expression consistent with 
the Kramers- Kronig relations can be used. Also the explicit expression of the reservoir field 
operators Y (x, oj) in terms of the bosonic operators / (x, oj) are not needed in the following. 

Using the commutation relations (||,f|) it can be shown fli6| , |20"f that the field operators 
A (x) and E (x) and the polarization field operators X (x) satisfy the well known canonical 
commutation relations 

A(x),E(x')]=-^-5(x~x'), (14) 
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X(x),E fx') 



0. 



(15) 



X(x),A (x') 



Using the Green function flTTD, Eq. may be written as 



A (x, uj) = (x, uj) + (x, w) 



(16) 



(17) 



where 



(x, = / dx' G A (x,x',u) f (x',u) = -ia^ \ dy e ik{x y) f (y , uj) 



and 



(x, uj) = I dx' Ga (x, x', cj) / (x', u;) = —ia^j— j dy e 



-ik(x-y) 



;is) 



(19) 



denote vector potential operator components describing wave propagation in the positive 
and negative x-direction respectively. From Eqs. (fL8|, I9|) we easily find 



d x A^ (x, uj) = ikA^ (x, uj) — iaJ—f (x, uj) 



(20) 



d x A(^) (x, uj) = —ikA^ (x, uj) + ia^j—f (x, uj) 



(21) 



These relations can be considered as quantum Langevin equations governing the spatial 
evolution of the vector potential operator components A^ (x, uj) and (x, uj). Here the 
operators / (x, uj) play the role of the Langevin force operators. Note that these noise terms 
are proportional to ^J~E\ which is responsible for the absorption. 

A similar decomposition as in Eq. (|T~7|) can be performed for the electric field operator 
components, where now the Green function Ge (x, x', uj) = iujGa (x, x', uj) has to be used. 
Finally we consider the time evolution of the operators in the Heisenberg picture. We 
immediately find from Eqs.(|l]^) 



id t f(x,uj)= f(x,uj),H L =uf(x,u) 



(22) 



and similar equations for the other field operators, which are connected with / (x,u). E.g. 
using Eq. ([i~8|) we have 



id t A^ (x, u) = A^)(x,uj),H l = U A(^)(x,u) . 



(23) 



It is seen from Eqs. (p2 , p3| ) that the operators show a harmonic time evolution in the case 
of a linear dielectric: 



f(x,uj,t + t) =e- wt f(x,w,t Q ) 



(24) 



A M (x, u, t + t) = e-**il H (s, w, t ) , (25) 

and u can be associated with the frequency of the harmonic time evolution of all the oper- 
ators with the argument uo. Thus, in case of a linear dielectric the decomposition in Eq.(fj|) 
can be considered as a Fourier decomposition of the operator of the vector potential in the 
Heisenberg picture. 

Concluding this section we state that the formulae given above enables one to study 
quantum properties of light fields propagating in a linear dielectric with absorption and 



dispersion |24| or passing through a multi-layer structure |23|. 



III. NONLINEAR CASE 

In this section we extend the linear model of a dielectric in order to include non-linear 
optical effects. This offers the possibility for studying the propagation of quantum light 
fields in nonlinear optical media including dispersion and absorption. In the Hopfield model 
@] the dielectric is described by a polarization field, the equations of motion for this field 
being linear. In reality the material system is not composed of harmonic oscillators and so 
this description is only an approximate one. In order to introduce into the Hopfield-model 
the really existing nonlinearities of the dielectric we add to the Lagrangian of the material 
system an arbitrary functional — $[X {x)\ of the polarization field X (x). This term does not 
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affect the definition of the canonical momenta conjugated to the field variables. Therefore 
the whole quantization scheme developed for the linear dielectric can be used also in the case 
of the nonlinear dielectric. The Hamiltonian describing the nonlinear dielectric interacting 
with the radiation field then reads as 



H = h {n L + n N ) = h J dx (H L (x) + H N (x) 



(26) 



where Hl (x) is the linear part of the Hamiltonian density which is defined in Eq.(|2|), and 
Hn (x) describes the intrinsic nonlinearity of the dielectric, 



H N (x) = $[X (x) 



(27) 



We emphasize that all the equations (||) - (|2ll) between the electromagnetic field opera- 
tors, the polarization field, and the bosonic creation and destruction operators f> (x, u) 
and / (x, to) are valid also in the case of the nonlinear dielectric interacting with the elec- 
tromagnetic field. These equations have been derived in the case of the linear dielectric 



by the diagonalization of the linear Hamiltonian ||T6| , |20[| , but without using the Heisenberg 
equations of motion for the corresponding operators. Therefore in the case of the nonlinear 



dielectric described by the Hamiltonian (|26| , |2"7|) , the same procedure can be applied. But 
with respect to the time development in the Heisenberg picture instead of Eqs.((2^ 
we find 



now 



id t f(x,uj) = f (x, uj) , H L + H N =u f(x,u)+ f{x,uj),H 



IN 



(28) 



id t A(^ (x, u) — (x,uj) ,Hl + Hn =wi(^) (x, uj) , H 



IN 



(29) 



These equations of motion, together with the quantum Langevin equations ( pO| , pl|) governing 
the spatial evolution of the vector potential operator components describe the space time 
behavior of the radiation field in a non-linear, dispersive and absorbing dielectric. These 
equations have to be supplemented by the relation (|Tt]) between the polarization field op- 
erator X (x, uj), the operator of the electric field E (x, ui) and the bosonic operator / (x, uj). 
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Using this equation we can express the non-linear part of the Hamiltonian completely in 
terms of E(x,u) and f(x,u) and their Hermitian conjugate operators. As a consequence 
Eqs. (|28|j29| ) include nonlinear optical effects in the time development of the radiation field. 
As can be seen from the Heisenberg equations of motion (|28|J29| ) due to the nonlinear part 
of the Hamiltonian the time development of the corresponding operators is much more com- 
plicated as in the linear case. Therefore the decomposition Eq.([|) now cannot be considered 
as a Fourier decomposition of the operator of the vector potential in the Heisenberg picture, 
because these components show a very complicated time behavior, and the parameter to 
cannot be associated with the frequency in the Fourier space. 

Finally we would like to note that the diagonalization procedure remains the same also in 
case of a more general nonlinear part (x) of the Hamiltonian, which may be a functional 
of the operators A(x), X (x), Y (x, uo) and its spacial derivatives (see e.g. f2l|), and the Eqs. 
(0) - (^1|) as well as the Heisenberg equations of motion fl28[ , |29"D may be used in order to 



investigate the radiation field inside the nonlinear dielectric. This offers the possibility to 
generalize the Hopfield model if necessary. 

IV. PULSE PROPAGATION EQUATION. 

In order to study the propagation of quantum pulses in nonlinear dispersive and ab- 



sorptive dielectrics we have to solve the nonlinear propagation equations in time (p8| 



together with the quantum Langevin equations ( pOjpTD governing the spatial evolution. In 



this section we shall consider these equations in more detail in the case of narrow-frequency 
band quantum pulses, and eventually for a Kerr medium we derive the quantum version of 
the soliton equation including also the effect of absorption. 

To begin with we divide the u -interval [O..00] into subintervals and decompose the field 
operators as follows: 

00 

A«(*) = £A (+) (*), (30) 

i=0 
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Aw/2 




(31) 



-Aw/2 



where 



u % = (i + 1/2) Au. 



(32) 



The analogous relations can be written for all the other operators which can be expressed 
by an integral over u. Now we assume that the light pulse under consideration has a small 
spectral bandwidth, so that only components of the field operator in one of the intervals, 
characterized by an index i , must be considered. It means that the bandwidth should 
not exceed the length of the corresponding sub interval Au. This also gives a criterium for 
choosing the value of Auj. In the following we restrict the parameter u to this interval, and 
we write 



(33) 



where Uq = Next we introduce slowly varying operators A/-*) {x,Q) (operators corre- 
sponding to the other direction can be introduced in a similar way) and /, (x, Q) according 



to 



A H) (x,u) = e l ^A^ ) (x,n), 



(34) 



(the index i will be suppressed here and in the following), where 



(35) 



and k m is a real number to be determined later. 



Using the Heisenberg equations of motion (|28], |29[) we find 




(36) 




(37) 
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where the notation 



Hn, O 



(38) 



has been introduced. We rewrite also the quantum Langevin equation ( pO]) in terms of the 
new slowly varying operators (Eqs. Q3H])): 

(39) 



d x A(^) (x, Vt) =i(k- k v ) (x, SI) - iaJ — / ( _^ (x, fi) . 



To be able to treat these equations (^,|^,^) in more detail we expand the wave number k 
as a function of to near uq into a Taylor series 



and in a similar way we have 



oo Si 171 



Si ^ W 



(40) 



E p 

m=0 



ml 



(41) 



where the coefficients k m = fc mr p m = p mr + zp m j are in general complex. Substi- 

tuting these expansions into Eq.(^) and employing Eqs.( |36| , |3~T|) , we arrive at 



\m=0 



ml 



d x A(->) (x, = i E -7 ( tdt + H *) ~ k f n ) 



m=0 



Integration over f2 gives 



9a;A|^) (x) = i 



Vm=0 



where A(^) (x) is given by the relations 



(42) 



(43) 



i(+) (a;) = i[+ ) ) (x)+i|t ) )(^), 



(+) 



iW ( x ) = (x) 



(+) 



(44) 



and 
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Auj/2 



R.S. = -ia ]T fa + H x N ) m J dQ (x, fi) 

m=0 m " A . . /n 



(45) 



-Aw/2 



Note that the various operators A and / should be labeled by the index io, describing the 
parameter range to which these operators belong (see e.g. Eq. Q50|) ) . The equations 
- |45|) describe the propagation of a quantum pulse with a narrow bandwidth spectrum in 
the dielectric with any order of dispersion, absorption and nonlinearity. In general we need 
similar equations for the fields running in the other direction because they are included in 
ii-N- It should be mentioned that the restriction to a pulse with narrow spectral bandwidth is 
not necessary, but in the general case we must consider for each u -interval such an equation 
as Eq.©. 

Let us consider a special case with small dispersion, where the absorption in the frequency 
interval around u io can be neglected (e* (cu) =lm(e (a;)) = 0, k mi — , k± r = l/v gr , where v gr 
is the group velocity, and k mr = for m > 2). Choosing k v = k 0r we arrive at the simplest 
case of a propagation equation for a quantum light pulse in a nonlinear dielectric, 



d x H d t 



-gr 



Km, (x) 



(46) 



Note that the assumption of zero absorption resulted in the vanishing of the term pro- 
portional to the operator /, , (x,Q), which may be considered as a Langevin fluctuation 
operator. Such operators have to be included in the Heisenberg equation of motion in order 
to describe correctly also losses. 



A. Second order dispersion approximation. 

Now we consider the propagation of a quantum pulse with narrow spectral bandwidth 
including to some extend dispersion and absorption. We neglect all terms proportional to 
Q m with m > 2. That means we take into account the dispersion only up to the second order. 
Moreover we neglect all terms in which the nonlinear Hamiltonian Ti^ and the Langevin 
fluctuation operators f , , (x, Cl) are coupled with higher order dispersion and absorption 

12 



effects (this means terms proportional to kmH^, m > 2 andp m jf, (x, Q),m> 1). Choosing 
k<p = k 0r , from Eq.(f43|) we arrive at the following propagation equation for a quantum light 
pulse 

(id x + ik oi + ihd t - l/2k 2 d tt + hUff) A[5) (x) 

Aw/2 

= « Po | dfl/^fofi), ( 47 ) 

-Aw/2 

which will be discussed in connection with a Kerr medium in the next subsection. 



B. Approximation for the propagation in a Kerr medium. 

In this section as an application of Eq.(fI7|) we discuss the pulse propagation in a Kerr 
medium. At first we specify the nonlinear part of the Hamiltonian H.n- The most simple 
extension of the linear dielectric to include nonlinearities is in case of a nonlinear dielectric 
with inversion symmetry the Ansatz 

*[*(*)] = ^jX*(x), (48) 

from which we find (see Eq . fl27|) ) 

H N = - J dx'X A (x f ) , (49) 



-oo 



where A is a constant characterizing the nonlinearity in a homogeneous dielectric. In consid- 
ering the propagation of a quantum pulse with narrow spectral bandwidth into the positive 
x-direction (see Eq. (|47|)) we may rewrite Hn approximately as 



Hn = j[ dx' X_[2) (x') (x') IJ^ (x') 1£) (x') , (50) 



where (remember Eq.fl33|)) 

Aw/2 

(x r ) = J d,tt [iu {e (u) - 1} 2 H) (x, fi) 



-Aw/2 



-2iac^Ei (u)f,s (x, f2) 
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(51) 



and the operators A(^) (x,ft) and (x, Q) are defined in Eq. fl3~4l) . Here we have also 

introduced a normal ordering of the operators and quadratic terms in X have been neglected. 

,(+) 



Straightforward calculations give for the commutator term 



of Eq. 



- J dx' G H (x, x') (x') X^ (x') ij^ (x') , (52) 



where 

Aw/2 

G H (x,x') = a 2 - I dn[iu{e*-l)-^ n -2c^U{Ax)) 
Pj n/2 \ 2ki\e\ J 

x exp [—k\ | Ax | + i (k T — k 0r ) Ax] , (53) 

U (Ax) denotes the unity step function, and Ax = x — x'. We note that the relation 
in general describes a nonlocal, nonlinear response of the matter, and Eq. ( fT7| ) is an 
integro-differential equation with respect to the space co-ordinate x. 

Now we consider the limiting case of weak absorption in the interval around ujq. Em- 
ploying the approximation 

Si < e r , (54) 

the term may be approximated by 



£i £i £i 1 



(55) 



2M|£| 2^I mv ^T7f|£| 2^Im(l+^)£ r V 
and the term containing the step function can be neglected, because it is proportional to 
As a further approximation we expand the integrand of the integral kernel in powers of Q 
and retain only the lowest order terms arriving at 



Alo/2 

Gh (x, x') = a 2 — iujq (e* — 1) — / dVt exp [iQ Ax fci r ] , 

Kr J 

' Aw/2 



2 C 0, . 1N 2 



= a 2 -iuo (e* - 1) — — — — sin [Aw k lr Ax/2] . (56) 
p v 1 k r Ax k lr L v 7 

Assuming that the operators X^ui) in Eq.(|52"D are slowly varying on a scale defined by the 
characteristic length [Auj /ci r ) _1 of Gh (x,x f ), we can replace the function Gh (x,x f ) by a 
delta-function, and the commutator term in Eq . (|52"D reads as 
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H N ,A^ (x) 



-iu (e - 1) X H) (x) X H) (a;) X H) (a;) . (57) 



Higher order approximations in Eq. (|53D give rise to nonlocal corrections in Eq. (0), which 
will be neglected here. In the same way as we have arrived at the Eq.([|7]), we neglect in 
Eq.flSTD the Langevin fluctuation operators /, , (x, Q), which result from the decomposition 
of (x) (see Eq. ([)]])), and we substitute k lr for k\, thus finding eventually 



' lr 



X Ah) (x) 4?) (s) (*) , (58) 



where 



^^(> |£ - 11 ) 4 ' (59) 

Substituting this in Eq. ([47]) we obtain finally the following equation for the propagation 
of a narrow-frequency band quantum pulse in a Kerr medium with dispersion and weak 
absorption 

id x + ik 0i + %k x d t - l/2k 2 d tt + % A[J } (x) (x)\ (a;) 

Auj/2 



= aJ- J dQf^ix,^. (60) 

-Aoj/2 

The absorption properties of the dielectric are described by the terms proportional to k mi 
and the Langevin fluctuation operator on the right hand side of this equation. The term 
proportional to \ describes the Kerr nonlinearity, and the dispersion effects are contained 
in the term proportional to kit- Of course the solution of this nonlinear operator equation 
containing also a Langevin fluctuation operator will be difficult, and numerical methods 
must be employed. 

It is clearly seen that in the case of pure real coefficients k\ and k% in Eq.(|60"D and vanishing 
E\ we obtain the operator version of the classical (see e.g. [0,0) Schrodinger equation similar 
to the one used in fll||l~4]| . We note that such an approximation is in general not consistent 



with the Kramers-Kronig relations because dispersion is always accompanied by absorption. 
Including the absorption brings an additional fluctuation term into the Eq.flBTiD which is 
important, in order that the canonical field commutation relations Eq.(0) are preserved. 
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V. SUMMARY 



On the basis of the quantization of the radiation field in a linear dielectric with dispersion 
and absorption as it has been described by Huttner and Barnett and Gruner and Welsch 
f20|| , we have developed a quantization scheme which includes besides linear dispersion and 
absorption also nonlinear optical effects. This has been achieved by introducing into the 
matter part of the microscopic Huttner-Barnett model a nonlinear interaction term depend- 
ing only on the polarization field operator which describes the dielectric. This guaranties, 
that the operator relations at equal times of the linear theory are also valid in the extended 
theory. Especially the equal time commutation relations of the free field case are preserved, 
and the space propagation equations for field components are the same as in the case of a 
linear dielectric. We have applied the theory to a narrow- frequency band quantum pulse 
propagating in a dielectric with a Kerr-like nonlinearity. Using a second order dispersion 
approximation and assuming weak absorption in the frequency band under consideration 
we have derived in the Heisenberg picture a nonlinear operator equation which governs the 
space-time development of the quantum pulse. It contains the Kerr nonlinearity, second 
order dispersion and an absorption term together with a corresponding fluctuation opera- 
tor. The resulting equation can be considered as a generalization of a nonlinear Schrodinger 
equation describing quantum solitons. 
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